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Abstract 

The linear stability of stalled accretion shocks is investigated in the 
context of core collapse of type II supernovae. We focus on a particular in- 
stability mechanism based on the coupling of acoustic perturbations with 
advected ones (vorticity, entropy). This advective-acoustic cycle takes 
place between the shock and the nascent neutron star. Both adiabatic 
and non-adiabatic processes may contribute to this coupling, but only 
adiabatic ones are considered in this first approach. The growth time of 
the adiabatic instability scales like the advection time, and is dominated 
by low degree modes 1=0,1,2. Non radial modes (1=1,2) found unstable 
by Blondin et al. (2003) can be related to this mechanism. 

1 Introduction 

Shocked accretion onto the surface of a compact star is known to be 
unstable in the context of magnetized white dwarfs, leading to shock os- 
cillations (from Langer, Chanmugam & Shaviv 1981, hereafter LCS81, to 
Saxton & Wu 2001). Houck & Chevalier (1992, hereafter HC92) made 
a linear stability analysis of shocked accretion onto a neutron star, and 
found an instability reminiscent of the instability found by LCS81. HC92 
showed specific cases where the cooling occurs mostly in a thin layer at the 
surface of the neutron star, while the flow is essentially adiabatic above 
it. The mechanism of the instability was described by LCS81 and sub- 
sequent authors as a kind of thermal instability: if the shock surface is 
moving outwards, the higher incident velocity in the frame of the shock 
produces a higher temperature blob, which pushes the shock further out if 
the increased cooling time exceeds the increased advection time. This cy- 
cle, however, resembles the unstable adiabatic cycle described by Foglizzo 
& Tagger (2000, hereafter FTOO), and Foglizzo (2001, 2002, hereafter 
F01,F02) in the context of shocked spherical accretion onto a black hole. 
In this case the acoustic feedback is purely adiabatic and is due to the ad- 
vection of vortical/entropic perturbations from the shock to the accretor. 
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A similar coupling must also take place if the accrotor is a neutron star. 
Is the instability found by HC92 due to a cooling process, in the spirit 
of LCS81, or is there a significant contribution of the adiabatic coupling 
between vortical/entropic and acoustic perturbations ? The recent work 
of Blondin, Mezzacappa & DeMarino (2003, hereafter BMD03) seems to 
support this second hypothesis. However, is the acoustic feedback identi- 
fied by BMD03 linear or due to turbulence ? Understanding the physical 
mechanism underlying this instability could prove useful in order to eval- 
uate its role when realistic non-adiabatic processes are taken into account. 



2 The perturbed adiabatic flow viewed 
as a forced oscillator 

In order to distinguish between adiabatic processes and non adiabatic 
ones, the flow structure between the shock Vah and the surface of the 
nascent neutron star is schematized as an adiabatic flow above a cooling 
layer rcooi. The present study focuses on the stability of the adiabatic part 
of the flow rcooi < r < Tsh. The entropy and vorticity equations can be 
integrated explicitly as in FOl. The Rankine-Hugoniot conditions at the 
shock impose that entropy and vorticity perturbations SS, 5w are simply 
related through: 



tJi X X ^ fn dSS dSS\ 
(5wr,5we,5w^) = 0, . ^ 

\ 7rt;sm6' 0(p jrv oo j 
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The differential system satisfied by perturbations is the same as in FOl 
(Eqs. (B18-B19)), only the functions A^, c describing the stationary flow 
are different. Pressure perturbations satisfy a differential equation with a 
source term due to entropy/vorticity perturbations (Eq. (4) of FOl) : 

f 9^ d ) 5p -A 5ps 

|a;3 + «ia;: + «o|- = _ — ■ (2) 

The source term 5ps is: 

^ ^ (3) 
p luj or \ \ iu> or J !\ ^ 

A = uj' + l{l + lf^ + 2ko^,. (4) 

Eq. (3) thus describes the local production of pressure perturbations due 
to advection of entropy and vorticity perturbations in a inhomogoneous 
flow. It characterizes the "excitator" , whereas the left hand side of Eq. (2) 
characterizes the "oscillator". According to Eq. (3), the local strength of 
the excitator depends on the value of the ratio uir/v. 

A f Sps\ _ 1 iivr dlogc^ 
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6S if 1 = 0, (6) 
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SS if « > 1. (7) 



In view of Eqs. (5) to (7) and the advective-acoustic cycles described in 
FOl, F02, two physical processes couple advected perturbations to the 
acoustic field: 

(i) The gradient of temperature characteristic of the entropic-acoustic 
cycle is essential for spherical perturbations 1 = (Eq. (6) and FTOO) and 
high frequency perturbations (Eq. (5) and FOl). 

(ii) Even in a isothermal flow (i.e. dc jdr = 0), non radial perturba- 
tions can excite acoustic waves in a vortical-acoustic cycle at low frequency 
(Eq. (7) and F02). 

The simple estimate in Eqs. (5) to (7) shows that the strength of the cx- 
citator is comparable for radial and non radial perturbations, and that its 
amplitude is highest near the lower boundary. 

An efficient coupling between the excitator and the oscillator requires not 
only a strong amplitude of the excitator, but also a good matching of 
their spatial lengthscales. The wavelength of the excitator (~ 2-kv/uj) is 
approximately a factor M smaller than the wavelength of the oscillator 
(~ 27r(c ± v)/u>). This contrasts with the simpler case of black hole ac- 
cretion, in which case the excitator and oscillator have comparable wave- 
lengths near the sonic radius {M. — > !)• The numerous oscillations of the 
excitator per acoustic wavelength should thus lead to a weak efBciency of 
the advective acoustic coupling in the inner regions where <C 1. The 
inner regions, however, are precisely the place where the amplitude of the 
excitator is highest, because the adiabatic gradients are strongest there. 
What is the net effect ? The choice of the lower boundary condition is 
crucial in answering this question. 



In order to separate the adiabatic effects from the non adiabatic ones, we 
choose to estimate the contribution of the adiabatic region by neglecting 
the acoustic feedback from the cooling layer as much as possible. The 
following assumptions are made at the lower boundary rcooi: 

(i) acoustic perturbations propagating downward are perfectly reflected 

out {to < c^ooi/r cool } 7 

(ii) entropy and vorticity perturbations are freely advected below rcooi 
in the cooling region, where their coupling to acoustic waves is ignored. 
Condition (ii) is equivalent to imposing that below rcooi, entropy and vor- 
ticity perturbations cease to be source terms of the acoustic equation. The 
source term in Eq. (2) is thus artificially damped by multiplying it by a 
smooth transition function <E>a. The transition is assumed to take place 
over a length A, comparable to the cooling length. This damping of the 
source term can either be viewed as an ad-hoc damping of the entropy 
perturbation itself, independent of its frequency, or as a smoothing of the 
flow gradient responsible for the coupling. 

The equations corresponding to these assumptions are obtained by match- 
ing the pressure perturbation 5p for r > rcooi with the homogeneous so- 
lution Spo associated to Eq. (2) for r < rcooi ■ 



3 Boundary condition at TcooI 
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4 Eigenmodes of shocked accretion 



The Rankinc-Hugoniot jump conditions arc used to compute the per- 
turbed quantities after the shock. These calculations are similar to those 
of Landau & Lifshitz (1987, Chap. 90) extended to the case of a non uni- 
form flow, or Nakayama (1992) extended to non radial perturbations. For 
a perturbed shock velocity Av in the strong shock limit. 
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The boundary value problem was solved numerically for I = 0,1. A broad 
range of unstable modes grow on a timescale comparable to a fraction of 
the advection timescale. The radial mode is always the most unstable, 
closely followed by the non radial mode 1 = 1. 



5 Evidence for the advective-acoustic cy- 
cle 

Following the same method as F02, the discrete spectrum obtained in the 
boundary value problem is checked by computing, in two steps, the effi- 
ciency Qadv of sound production by the advection of an cntropy/vorticity 
perturbation (without a shock), and the efficiency Qsh of entropy /vorticity 
production by an outgoing acoustic wave reaching a shock. Qadv is defined 
by: 

The efficiency Qsh is obtained through a WKB approximation at high 
frequency: 

I6sh|= \ ^ ■ (13) 

Ml 1 + 

The global efficiency Q = | Qadv Qsh | of the advective-acoustic cycle leads 
to a first estimate of the growth rate at high frequency: 

Wi=— log Q, (14) 

Ttot 

where rtot is the total duration of the cycle (advection -|- acoustic). The 
exact resolution of the discrete eigenfrequencies was successfuly compared 
to the continuous WKB estimate (14). 
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6 Comparison with BMD03 



Our stability analysis in an adiabatic flow should apply directly to the 
numerical simulations of BMD03. The authors recognized the existence 
of an advective-acoustic cycle similar to the one occuring in the vortical- 
acoustic instability of F02. Our linear study seems to agree qualitatively 
with their results concerning the mode / = 1. Nevertheless, the mode 
Z = is stable in their simulations, whereas it is the most unstable one 
according to our calculations, as well as in the work of HC92. This im- 
portant difference may come from the "leaky" lower boundary condition 
of BMD03, which is different from ours even in the linear regime. 

7 Conclusion 

The temperature and velocity gradients in the subsonic flow between the 
shock and the accretor is responsible for an efficient linear coupling be- 
tween entropy /vorticity and acoustic perturbations. Most of the sound 
comes from the region close to the lower boundary of the adiabatic re- 
gion, despite the fact that the wavelength of advected perturbations is 
much smaller than acoustic ones. The acoustic waves reaching the shock 
produce new entropy /vorticity perturbations, in an unstable cycle. The 
growth time is comparable to the advection time. The most unstable 
modes correspond to Z = 0, 1 perturbations. The identification of the 
advective-acoustic cycle as the mechanism responsible for the instability 
was checked numerically through the calculation of Qadv and Qsh- Since 
the region contributing most efficiently to the instability is the vicinity of 
the lower boundary Tcooi, the role of cooling processes is at least crucial 
in determining the growth rate of the instability. The principle of a cycle 
between propagating and advected perturbations could be useful in inter- 
preting the stability analysis of non adiabatic flows. Whether non adia- 
batic processes are partially stabilizing or even more destabilizing remains 
to be determined. The difficulties of 1-D numerical models in reaching an 
explosion (e.g. Buras et al. 2003) suggests that cooling processes are 
indeed important enough to significantly stabilize the entropic-acoustic 
cycle at work for radial perturbations. 
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